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Recently Kaloper, Kleban and Martin reexamined the McVittie solution and argued, contrary to 
a very widely held belief, that the solution contains a black hole in an expanding universe. Here 
we corroborate their main conclusion but go on to examine, in some detail, a specific solution that 
asymptotes to the ACDM cosmology. We show that part of the boundary of the solution contains 
the inner bifurcation two - sphere of the Schwarzschild - de Sitter spacetime and so both the black 
and white hole horizons together form a partial boundary of this McVittie solution. We go on to 
show that the null and weak energy conditions are satisfied and that the dominant energy condition 
is satisfied almost everywhere in the solution. The solution is understood here by way of a systematic 
construction of a conformal diagram based on detailed numerical integrations of the null geodesic 
equations. We find that the McVittie solution admits a degenerate limit in which the bifurcation 
two - sphere disappears. For solutions with zero cosmological constant, we find no evidence for the 
development of a weak null singularity. Rather, we find that in this case there is either a black hole 
to the future of an initial singularity or a white hole to its past. 
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I. INTRODUCTION 



A simple, but painful truth is the fact that it is far easier to find an exact solution to Einstein's equations than it is 
to understand it. A fine example of this is given by McVittie's inhomogeneous cosmological solution the meaning 
of which has been debated since 1933. In retrospect, this effort has to be considered an utterly remarkable step into 
an area of research which is, to this day, still in its infancy Q. The McVittie solution has been the subject of a large 
number of investigations (we point to the recent thesis by Martin Q and the review given in Q ) , and generalizations 
0, but only recently did Kaloper, Kleban and Martin ^8] (henceforth KKM) explain the misinformation which has 
developed around this solution. In this paper we corroborate the main conclusion in KKM but also penetrate more 
deeply into an understanding of a specific solution. By way of the specification of a characteristic function for the 
solution, we exhibit a specific solution that asymptotes to the standard ACDM universe. Remarkably, part of the 
inner boundary of this solution contains the inner bifurcation two - sphere of the non - degenerate Schwarzschild - de 
Sitter spacetime. This tells us that both the black and white hole horizons of the extended Schwarzschild - de Sitter 
spacetime form part of the boundary of this McVittie solution, a possibility not envisioned in the KKM analysis. The 
specification of a characteristic function provides sufficient detail to allow for an examination of energy conditions 
and a systematic construction of the conformal diagram based on detailed numerical integrations of the null geodesic 
equations. The present work suggests that the very notion of an inhomogeneity in cosmology may well go beyond 
the concept of inhomogeneity in elementary physical variables. Finally, a minor point in the KKM analysis was the 
suggestion that in the case of a zero cosmological constant, the would - be horizon forms a weak null singularity. Here 
we find no evidence for this behavioural 



II. THE SOLUTION 



Overview 



The particular solution we are concerned with here is the simplest of the McVittie classic and this can be written 
in the form (e.g. tlfl]) 



' dt^ + e/3(*)(i + m/2u)\dr^ + r^dnl) (1) 

1 + m/2u J 



where u = re^^'^, m is a positive constant and dill is the metric of a unit 2-sphere. Clearly, for m = we obtain 
a spatially fiat Robertson - Walker metric, and for constant (3, we have the Schwarzschild metric (here in isotropic 
coordinates). These observations do not constitute an understanding of the metric Indeed, it is a trivial exercise 
to construct distinct spacetimes with these two fundamental features. Perhaps, the enduring interest in the McVittie 
solution derives from the observation that if we take the coordinates of ([T]) as comoving then we obtain a perfect fluid 
with energy density p and isotropic pressure p given by (e.g. [lOj ) 



where ■ = d/dt and 



4'^ ' ^ 4^^ ^l-2m/R{t,r) ^ ' 



R = u{l + m/2uf. (3) 



The uniform nature of the energy density and non-uniform nature of the pressure is often brought forward as a reason 
to consider this solution unphysical. However, even in the static Schwarzschild interior solution, such conclusions can 
be considered hasty [l3|- Our purpose here is not to argue, ab initio, for the physicality of the McVittie solution, 
but rather our purpose is to exhibit in detail the rather remarkable, and heretofore unrecognized, geometric structure 
that this solution presents. We comment on the idealization that the solution represents only once this structure is 
developed. 



^ Whereas the analysis given here parallels, in some respects, that given in KKM, it also differs in a number of important aspects. We 
present a full discussion and point out at various stages agreement and disagreement with the KKM analysis. 

^ Recently, study of the McVittie solution has been denigrated 3 with the view, in part, that this class of solutions is but a simple subset 
of a larger class of known solutions. The results presented in this paper argue for the contrary view. Even a simple looking metric can, 
when properly studied, yield a rich geometric structure. 
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Under the coordinate transformation defined by ([3]), the metric ([T]) becomes (e.g. (Toj ) 

2H{t)R , 
^1 - 2m/R 1 - 2m/i? 



where 

f EEl-2m/R- H^R^ (5) 

and is the Hubble function, given hyH — $/2^d/a where a{t) is the usual scale factor for to = 0. The form (|4]) 
is the basis for much of the analysis in KKM. Here we use the form ^ and additional transformations suitable for 
numerical integrations. Let us note that the effective gravitational mass |13|] associated with ^ is not to, but rather 
M, given by 

M{t,R) =m+^H^R^. (6) 
B. The Function H 

Here we are not interested in arbitrary functions H , but only those that reflect, in a general way, the background 
cosmological model as it is currently understood. In particular, we take < for finite t, 

lim H = Ho>0, \im H ^ hm H ^0, \im H = oo, (7) 

and 

lim /3 = -00. (8) 
These general properties are in fact crucial to the present analysis. 



As recognized in KKM, i = is not, in general, part of the spacetime. To see this here, from the definition of u 
and ([8]), we find limt_>o u = for all finite r. From the definition ([3]) then we have 

^ fO ifTO = 

lim R^ i , (9) 

t^o 1 cx) if TO 7^ 

and so f = is not, for to 7^ 0, part of the spacetime (and no limit to — exists). This result depends on ([5]) and 
relaxation of ([8]) gives rise to other more involved possibilities not discussed here. 



C. Scalar Singularities 

Let us consider the singularities of as revealed by scalars polynomial in the Riemann tensor. It turns out that 
we need only report the Ricci scalar TZ, which we find is given by 

7^ = 12H^ + , (10) 

^l-2m/i? 

in agreement with KKM, since all other invariants, derived from (partial) derivatives of the metric tensor no higher 
than 2, add no new information. For < t < 00, there is a singularity at i? = 2to {u = to/2), which, as is clear 
from ([1]), is spacelike in agreement with KKM. The apparent singularity at t = 0, over the range 2m < i? < cx) is, as 
explained above, not part of the spacetime (|1]) due to ([5]). 
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D. Asymptotics 



With conditions ([7]) let us also consider the asymptotic limit t — > oo, and in particular the roots to /o = 1 — 2m/ R — 
HqR^ — 0. There are three cases: 21m^Hf^ > 1, for which there are no positive roots, 27m^7?g = 1, for which there 
is one (coincident) positive root R = 3m, and 27m^HQ < 1 for which there are two distinct positive roots that satisfy 
< 2m < i?_ < 3m < It is only the last case which is of central interest here. As explained in detail below, we 
are interested in solutions that asymptote to de Sitter space for R and t ^ oo and become Schwarzschild - de Sitter 
space for R -> _R_ and t — > oo. This requires (e.g. [13]) 



A/3 



where A is the cosmological constant. For 27m^i7g <C 1 we note that i?_ ~ 2m and R^ 
(e.g. [l^). The solution discussed here is essentially homogeneous at this length scale. 



(11) 

I/Hq, the Hubble length 



E. The Locus / = 

Tangents to the locus / = 0, where / is defined by are null for R = 2H'^R^, that is, H = 2H{3m - R)/R'^, 
and timelike for R < 2H'^R^. The tangents are spacelike for R > 2H'^R^ and for 7? < 2H'^R^ along the branch 
dR < 0,dt > 0. The essential features are summarized schematically in Figure [TJ 
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FIG. 1. The R — t plane and the locus / = 0. The locus is timelike above and to the right of n where H — 2H{3m — R)/R^ . 
The locus is spacelike below n. R — 3m is also shown. It passes through / = at o where t = T such that 27m^ H{T)^ — 1. 
R = 3M is spacelike to the left the locus and timelike to the right. The trajectories R± are spacelike. The asymptotic points 
are defined by a : — > oo, i? = 2m), b: (t — >■ cxa, i? = _R_) and c: {t oo,R = R+). These are explained below. 



III. NULL GEODESICS - QUALITATIVE 



A. Outgoing and Ingoing Geodesies 



Let us first examine, qualitatively, general properties of the radial null geodesies {r]) of These must satisfy 



dR 
It 



y/l-2m/R[HR± Vl - 2m/i? 



(12) 



in agreement with KKM. Now, whereas for the "+" ("outgoing") branch clearly dR/dt > 0, the "-" ("ingoing") branch 
requires further examination. From ([S]) and it follows that for the "-" branch dR/dt > for / < 0, dR/dt < 
for / > 0, and dR/dt = for / = 0. In particular, note that for / < 0, dR/dt > along both branches of ([12]). It 
is appropriate here to note from Q that tangents to surfaces of constant finite t are spacelike for R > 2m (and so 
for finite t we set the future orientation dt/dX\rf > for affine A increasing to the future) and tangents to surfaces of 
constant R are spacelike for / < 0, null for / = and timelike for / > 0. 
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B. Expansions 

Letting k^. signify the 4-tangents to the radial nuh geodesies, the associated expansions 9± = Vafc± follow as 



0± = 




HR± 




(13) 



where evaluation along T] is now understood. This is in agreement with KKM up to the last term which is missing 
in the KKM analysis. The change in sign of 6*- at / = led KKM to refer to / = as an "apparent horizon". 
Normally, this description would be reserved for a change in sign of 0+ since a change in sign of 0_ does not hide 
events below / = from distant observers. Now to see the importance of the last term in consider t oo along 
ingoing geodesies. Clearly all but the last term — > 0. However, as we show below, f — oo at finite A and so the last 
term in (jl3p diverges and as a result, 0_ becomes indeterminate in these coordinates. The removal of this ambiguity 
is discussed below. 



C. Infinity 



Let us examine the "outer" boundary of (|1]): {R 
equations can be written in the form 



oo, t 



3) = Z+ (spacelike [15|). The radial null geodesic 



dt 
dX 



< 



(14) 



where the inequality holds for finite t. From it follows that R ^ 00 for finite A and so both branches are future 
null geodesically complete 0. This is no surprise since I"*" is indistinguishable from 1+ for de Sitter space (given 
(1111) '). The "inner" boundary of Q requires a much more detailed analysis. This is explained below. Here we simply 
introduce the notation: (i? = R^,t = 00) = % (null), and note that 



and (i?_ < R < 00, t = 00) — 



/|„ =(i?-2A/)|„ =0, 



IV. A SPECIFIC FORM FOR H 



(15) 



In order to examine the solutions to (112]) we must first specify a specific background function H, which satisfies the 
conditions (O and (|S]), since (fT^ appears to have no analytic solution for general H. Here we choose 



H = 



Ha smh{3Hat) 



cosh(3iJoi) 

so as to reflect an asymptotic ACDM universe]! 



= iJocoth(— — ) 



(16) 



V. ENERGY CONDITIONS 

Note that we apply the classical energy conditions to ([T]) by way of the Einstein equations (without an explicit 
cosmological constant). As is clear from ([2]), with our general conditions on 77, p > and p + p > and so the null 
and weak energy conditions are satisfied. For the dominant energy condition, —p < p < p and so from ([2|) 
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— < 7/3'- 
2m /R 4 



(17) 



^ The same conclusion has been obtained by Brien Nolan (private communication). 

* This choice does not limit all of what follows. For example, the completeness/incompleteness arguments given in Appendix A are 
unaffected by this choice. 
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Whereas the left hand inequahty is always satisfied, given our general conditions on _ff , the right hand side can be 
given in the form 



where, from (fTBj) , 



R{t) > ■ip{t)m 
2(cosh(3ffoi) + 1)^ 



(18) 



(19) 



cosh{3Hot){cosh{3Hot) + 2) 

Now whereas the explicit value of ip at any t depends on Hq, the general form of ip does not. In particular, limt_i.o V' = 
8/3 and limj^cio "0 = 2. The function ip is shown in Figure [2] We conclude that the dominant energy condition is 
satisfied almost everywherc0. 



8/3 



io 



12 



14 



16 
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FIG. 2. The function tp given by (|19|l . The curve has been constructed with Hg — 2.3 10 ^. T is defined hy t = lO'^ where 
[t] = s. The current epoch is shown as a dot. 



Finally, the strong energy condition requires p + 3p > 0. From ^ it follows that this condition requires 



From (|16p we find 



and so we arrive at 



2m ij2 
R(t)- 



^ _ (cosh(3i7oi) + 4) (cosh(3i/ot) - 2) 

^ - (cosh(3i/ot) + 4)2 



< 


t < 


to 


= 


t = 


to 


> 


t > 


to 



where 



to 



ln(2 + V3) 



2Hn 



(20) 



(21) 



(22) 



(23) 



and so the strong energy condition is satisfied everywhere for t < tQ. However, for t > the strong energy condition 
is satisfied only for 



where, clearly, 



5 = 



R{t) < 6(t)m 
2(cosh(3iJoi) +4)^ 



(24) 
(25) 



(cosh(3i?o0 + 4)(cosh(3Hoi) - 2) ' 

Since, as it is easy to show, i5 drops rapidly from oo at to to (limt^^oo S —)2, we conclude that the strong energy 
condition is eventually satisfied almost nowhere, as expected. 



^ We have no explanation for the curious position of the current epoch in Figure [2] 
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VI. INTEGRATION OF THE NULL GEODESICS 



A. Integration of (|12p in the R~ t plane 



We now examine numerical solutions to the null geodesic equations (fT2|) . subject to the choice (|16|) . First, consider 
the outgoing geodesies given by "+" in (IT^ . The integrations are shown in Figure [31 These geodesies are monotone 
in the R — t plane. In contrast, the ingoing geodesies ("-" in ([T2|)) show considerably more structure. This is shown 
in Figure ID We find 5 distinct types of evolution. Moving from the bottom right to the upper left in Figure |3] 
we find: (i) Geodesies which asymptote monotonically to H (again, defined by (i? = = oo)) without crossing 

the locus / = 0, (ii) A last geodesic that asymptotes monotonically to H without crossing / = (call it 7^2), (in) 
Geodesies which reach a maximum R < i?+ at / = and then asymptote monotonically to H, (iv) A last geodesic 
that asymptotes monotonically to R+ and terminates in i+(call it rji), and (v) Geodesies that cross i?_ and R+ and 
monotonically evolve to X+. Now the cases (i) and the limit (ii), are absent in the KKM analysis, but are central to 
our examination of the inner boundary of Since the form (|12p is particularly sensitive to error for these cases, we 
now introduce new coordinates to demonstrate, in more detail, that there do indeed exist ingoing null geodesies that 
reach H without crossing f — 0. 



2 m 















I / y y y 





FIG. 3. The outgoing solutions of the null geodesic equation (|12|l under the condition (|16p . The values taken are Ho = 1/3 and 
m = {1/{HoIq)) — 958041/2000000 ~ 0.479 (where lo is defined below). These values are of no consequence as we are simply 
interested in the qualitative behaviour of the solutions to p2p . The locus / = is also shown. The values of the roots R± are 
~ 2.29 and R- ~ 1.11. The cone shows the "leg" of the null cone under consideration. 



2 m 





FIG. 4. The ingoing solutions of the null geodesic equation (|12p under the same conditions as Figure |31 
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B. Integration of (|12p in the z — I plane 

We recast the problem as follows: Let 



z= 1 1 



2m 



and let 



We observe the ranges 



and 



Hm 



(i? = 2m) < z < 1 (R^oo) 



(26) 



(27) 



(28) 



(t -> 0) <l <lo = 



1 



Horn 



[t oo). 



(29) 



With the definitions (PB)) and (P7)) . equation ([T^ takes the form (using the asymptotic ACDM model as before) 

dz _ /I - z2 z{l- z2)2 



I 



/2 



3(Zg-P);- 



(30) 



for the ingoing case. At first sight this might not appear to simplify things, but it does. First note that we need only 
specify Iq. For the case of central interest {27w?Hq < 1) we have 



^0 > sVs. 



As regards initial conditions, we have 



which is regular over the range 



Moreover, we have 



dz 
dl 



/2 



z=Q,l=h 3^1(^0 '1) 



<li < lo- 



>0 /<0,?<^ 



(1-^2 



<0 f>0,l>^^. 



(31) 



(32) 



(33) 



(34) 



Numerical integrations are shown in Figure [5] Our conclusion is that null geodesies that reach H without crossing 
/ = are a fundamental feature of this solution. 



VII. GLOBAL STRUCTURE OF THE SPACETIME 



A. Construction of the Conformal Diagram 

The conformal representation of a spacetime (Penrose - Carter diagram [17]) is, of course, not unique in detail. 
However, all conformal representations must show the global structure of the spacetimeQ Here we construct the 



This useful definition for z was pointed out to us by Brien Nolan (private communication). 

All cases considered here are time-symmetric in the sense that the diagrams can be flipped upside down. Moreover, all diagrams can be 
rotated, interchanging the left and right hand sides. 
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0.2 0.4 0.6 0.8 1 

/// 



FIG. 5. Integration of {SO]) for lo = 2000000/319347 ~ 6.263 so that z- = 37/100 and z+ = -37/200 + 35893^''V200 ~ 0.76 
over the range < / < Zq. The limiting geodesies Tji and Tj2 are shown as is the locus / = (which connects H with i'^). The 
dashed curves connecting z — and H and i'^ with z = 1 indicate d?z/dl^ = 0. 

Penrose - Carter diagram in the following way: The construction of the diagram starts by solving the null geodesic 
equations numerically. This gives a general understanding of the global behavior of the spacetime. In the present 
case, all of the outgoing geodesies and all of the ingoing geodesies originate from the singularity at i? = 2m. These 
geodesies intersect with R — 2m at some finite value of t > 0. That is, any point in the spacetime {t > 0, R > 2m) 
can be connected to the past boundary {R — 2m, t > 0) by a unique null geodesic from each branch. We represent the 
boundary R — 2m as a horizontal line in a Cartesian plane (y = 0, — 1 < x < 1), setting the right end of i? = 2m at 
t — 0, and the left end at t = oo. To represent the interval < t < oo, from a; = —1 to 1, we use the transformation 
function 

t = A(l-x)^tan(^(l-a;)^), (35) 

where A and B are adjustable constants. A represents the value of t at the center of the line x — and B adjusts 
the position of b. These constants have no physical significance and where adjusted only for visual presentation (we 
chose A ^ 0.39 and B ~ 0.75). Once this line is adopted [3, any point in the spacetime can be projected onto the 
conformal diagram by solving for both null geodesies that pass through any point, and numerically solving for the 
value of t at which these two curves reach R = 2m; say ti, and t2- This is illustrated in Figure |6l Once these two 
values of t are found, we can represent them in the conformal diagram by solving for x in ([35]). This gives us xi and 
X2 ■ Finally, we find the conformal representation of the original point by finding the intersection of the same two null 
geodesies, but presented in the conformal diagram as y = —x + xi for one geodesic, and y — x ~ X2 for the other. 
This is illustrated in Figure [T] By our choice ([35]). all diagrams are strongly compactificd in the region R > . 



B. Null geodesies 

Under the procedure described above. Figure |3] is mapped into the right hand side of Figure |8] and Figure |4] is 
mapped into the left hand side of Figure [5] 



C. Surfaces of constant R and t 

Using the same procedure, at the left in Figure [HI we show surfaces of constant t and at the right in Figure [S] we 
show surfaces of constant R. The surfaces of constant t are spacelike for all finite t. The surfaces of constant R are 
spacelike for / < and timelike for / > 0. 
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2 m 



FIG. 6. Locating ti and t2 for any event via null geodesies. The locus / = is shown. Note that t increases to the right. 




X X 
2 1 

FIG. 7. Gonformal representation of the event {t, R). The locus / = is shown. Note that t increases to the left. 

D. The fluid streamlines 

The streamlines r = constant can be written out explicitly in 2; — Z coordinates in the form 

.(O^.a.h(ilog((|)(|^))) (36) 
where ^(^1) = 0. Transforming to the R ~ t plane, these streamlines take the form 

«W = 2™c«h(il„g(f|tf)) t") 

where R{ti) = 2m. These streamlines are shown in Figure [TUl Note that the streamlines do not cross H. This fact 
can be considered the source of relation ([T5|) and is the central part of the original construction [sj . 

E. The boundary and ingoing geodesies 

In Figure 1111 we summarize the boundary of @ and classify ingoing null geodesies which do not terminate on . 
We are now in a position do discuss the affine completeness/incompleteness of these geodesies. As is discussed in detail 
in Appendix A, all these geodesies are incomplete except r]i which we find to be complete. This incompleteness is the 
central point in KKM, but their analysis revealed only the section (&, i^) of Ti. Here we observe that r^2 terminates 
at & on "H and we observe very special properties associated with b as explained in Appendix B: 6 is characterized by 
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FIG. 8. At right, the conformal representation of the outgoing null geodesies as given in Figure [3] The locus / = is shown. 
The bottom horizontal line represents the singularity R — 2m. At left the conformal representation of the ingoing null geodesies 
as given in Figure [4] 

vanishing expansion for both the ingoing and outgoing radial null geodesies. This is the hallmark of a bifurcation two 
- sphere. This bifurcation two - sphere divides T-l into two sections: a "black hole" horizon to the future of b and a 
"white hole" horizon to the past of b. This is explored in the completion to the spacetime given below. 



F. A completion 

One possible extension of the McVittie spacetime, which is null geodesically complete, is shown in Figure fT2l where 
we have identified the inner bifurcation two - sphere of the Schwarzschild - de Sitter spacetime with that of the 
McVittie spacetime [1^. This enlarged spacetime now has a center of symmetry, the center of the Schwarzschild - 
de Sitter black hole at i? = 0. Ingoing null geodesies below r/i but above ?]2 terminate at the singularity i? = in 
the Schwarzschild - de Sitter black hole. Whereas 7^2 joins onto the "left" black hole horizon of the Schwarzschild - 
de Sitter spacetime, all "ingoing" geodesies below r]2 have R monotonically increasing, pass through the white hole 
horizon and terminate at i? = oo, that is, Z+ of the Schwarzschild - de Sitter spacetime. Whereas the degenerate 
case 27m^HQ = 1 is outside the cases of interest from a physical point of view, it offers a very instructive limit from 
a mathematical point of view since T]i and ?]2 then coincide and b disappears. This is examined in Appendix C. 



VIII. DISCUSSION 

We have performed a detailed study of a particular McVittie solution, by way of the specification of a characteristic 
function, that asymptotes to the standard ACDM cosmology and that contains an inner boundary that is a slice of 
the extended Schwarzschild - de Sitter spacetime. We have found that this inner boundary contains a bifurcation 
two - sphere where the expansion of both the ingoing and outgoing radial null geodesies vanishes. To the future of 
this bifurcation on this inner boundary we have found a black hole horizon at finite affine distance and therefore we 
corroborate the main conclusion in the work of Kaloper, Kleban and Martin. In addition, however, we have found a 
white hole horizon to the past of the bifurcation which is also at finite affine distance. In the degenerate limit of this 
particular solution the bifurcation two - sphere and black hole horizon disappear leaving only the white hole horizon, 
also at finite affine distance. We have shown that the null and weak energy conditions are satisfied and that the 
dominant energy condition is satisfied almost everywhere. The global structure of the solution has been constructed 
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FIG. 9. At left the conformal representation of surfaces of constant t. The locus / = is also shown. Moving from the bottom 
right to the upper left we have: to < t < T (where 27m^H{T)'^ = 1), t = T, t ^ ti > T and t2 > ti. The boundary has t oo. 
At right the conformal representation of surfaces of constant R. The locus / = is also shown. Moving from the bottom left 
to the upper right we have: 2m < Ro < R < R-, R = R-, R- < R < 3m, R = 3m, R = R+, Ri > R+. The right boundary 
has i? — >■ oo and the left boundary has R = i?_ . 




FIG. 10. The fluid streamlines r — constant > in the conformal diagram. The streamlines are timelike and r increases to the 
right. Note that the streamlines do not cross "H which is the limit r — ^ 0. 



systematically based on detailed numerical integrations of the null geodesic equations. In the case Hq = 0, the work 
of Kaloper, Kleban and Martin suggested that the would - be horizon forms a weak null singularity. In Appendix D 
we argue that this is not the case. Moreover, we argue that this case is rather less interesting than Hq > since the 
solutions can have a black hole horizon to the future of the singularity at i? = 2m or a white hole horizon to the past 
of i? = 2m. 

The present analysis relies on the (standard) condition ([5]). Relaxing this, it is clear that we can maintain the 
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FIG. 11. The conformal representation of the boundary to the spacetime (|4]) given (|16p . The boundary is defined as follows: The 
singularity R = 2m, 1+ = (i? = oo,t = oo), i+ = {R- < R < oo,t = oo), K- {R = R-,t ^ oo), ~i° = (2m < J? < -> oo) 

and +i° = (2m <R<oo,t>0). 




FIG. 12. An extension of the McVittie spacetime which is null geodesically complete. The inner bifurcate two - sphere of the 
Schwarzschild - de Sitter spacetime (6) is identified with that of the McVittie spacetime. 



conditions ^ but also violate ([9]). This means that in general the initial singularity will consist of two parts: the 
"pressure" singularity at i? = 2m, and a generalized "big bang" singularity a.t t = 0. In a sense this shows how 
"delicate" the solution is. The fact that the McVittie solution cannot represent a physically realistic inhomogeneity 
is, we think, best shown by a glance at Figure [10] and equation (fT5)) . By construction, no fluid streamlines can cross H 
and so the black and white hole horizons are present here by way of mathematical extensions, not physical processes. 
There is no easy fix for this situation, within the context of McVittie's approach, since a routine calculation shows 
that ([T|) is a perfect fluid if and only if dm/dt = 0. Nonetheless, the present analysis shows that a rather routine 
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looking spacetime, like can in fact harbor a rather exotic interior. Moreover, the present analysis suggests that the 
very notion of an inhomogeneity in cosmology may go beyond the concept of inhomogeneity in elementary physical 
variables. 

What we have done here can be expanded in a number of ways. First, of course, one could relax condition (|16p 
and consider a wider class of possibilities. In all cases one would find that if the vacuum boundary ("H) contains a 
bifurcate two - sphere, then this bifurcate two - sphere is also part of the McVittie solution itself. This geometric 
behaviour can be traced to McVittie's no - fiux condition which preserves the integrity of H. One can reasonably 
expect that the integrity of H is destroyed by any flux through it. To conclude, we believe that the McVittie solution 
is an instructive idealization very much like the Kruskal - Szekeres extension. 



ACKNOWLEDGMENTS 



This work was supported in part by a grant from the Natural Sciences and Engineering Research Council of Canada 
(to KL) and a Carmichael Fellowship (to MA) . We have several people to thank for discussions which have improved 
our paper. In the early stages we benefited from discussions with David Garfinkle and Matt Visser. As the project 
developed, we had extensive discussions with Damien Martin, Matthew Kleban and Brien Nolan. Of course, we do 
not mean to suggest that these people necessarily agree or disagree with what we have concluded here. 



[*] Electronic Address: lake@astro.queensu.ca 
[**] Electronic Address: majd@astro.queensu.ca 
[3] G. C. McVittie, Mon. Not. R. Astron. Soc. 93, 325 (1933). See also Ap. J. 143, 682 (1966) and General Relativity and 

Cosmology (University of Illinois Press 1962) 
[4] For discussions of a black holes embedded in an expanding universe see A. Krasiriski and C. Hellaby, Phys. Rev. D 69, 
043502 (2004) ;arXiv:gr-qc/0309119v2], P. Jacewicz and A. Krasiriski [arXiv:gr-qc/1101.5307vl] and references therein. For 
a recent discussion of inhomogeneous cosmologies see K. Bolejko, A. Krasiriski, C. Hellaby and M-N Celerier, Structures 
in the Universe by Exact Methods (Cambridge University Press, Cambridge, 2009) and references therein. 
[5] Damien Martin, Ph. D. thesis. University of California, Davis (2011) 
[6] M. Carrera and D. Giulini, Rev. Mod. Phys. 82, 169 (2010) [arXiv:gr-qc/0810.2712v2] 
[7] M. Carrera and D. Giulini, Phys. Rev. D 81, 043521 (2010) [arXiv:gr-qc/0908.3101] 
[8] N. Kaloper, M. Kleban and D. Martin, Phys. Rev. D 81, 104044 (2010) [arXiv:hep-th/1003.4777v3] 
[9] K, Bolejko, M-N Celerier and A. Krasiriski, Classical Quantum Gravity (in press) (2011) [arXiv:astro-ph/1102.1449] 
[10] B. C. Nolan, Classical Quantum Gravity 16, 1227 (1999). 

[11] We use geometrical units and a signature of +2 throughout. Functional dependence is designated only on the first appear- 
ance of a function. 

[12] C. W. Misner, K. S. Thorne and J. A. Wheeler, Gravitation, (W. H. Freeman, San Francisco, 1973) 
[13] M = TZg^"^ /2 where TZ is the Riemann tensor. 

[14] See, for example, K. Lake and R. Roeder, Phys. Rev. D 15, 3513 (1977). 

[15] R. Penrose in Battelle Rencontres, edited by C. M. DeWitt and J. A. Wheeler (W. A Benjamin, New York, 1968) 
[16] M. J. Jaklitsch, C. Hellaby and D. R. Matravers, GRG 21, 941 (1989). 

[17] J. Griffiths and J. Podolsky, Exact Space-Times in Einstein's General Relativity (Cambridge University Press, Cambridge, 
2009) 

[18] The choice of the shape of the singularity R = 2m as a straight horizontal line is arbitrary, as is the choice of the time 
transformation function (|35|) . Any spacelike curve (i.e. — 1 < slope < 1 ) with a finite length, and any function, that is 
monotonically increasing, and spans < t < oo in a finite interval of the new variable, could be used to construct the 
conformal diagram. 

[19] See, for example, E. Poisson, A Relativist's Toolkit: The Mathematics of Black-Hole Mechanics (Cambridge University 

Press, Cambridge, 2004) 
[20] K. Lake, Classical Quantum Gravity 27, 097001 (2010) [arXiv:gr-qc/1002.3600] 

[21] This is defined by £(2;)e^'^) = x. See, for example, R. M. Corless, G. H. Gonnet, D. E. G. Hare, D. J. Jeffrey, and D. E. 

Knuth, Advances in Computational Mathematics 5, 329 (1996). 
[22] With conditions Q the metric tensor and all derivatives to order 2 are continuous accross T-i and so matching accross T-i 

is complete under any formulation of matching conditions. 
[23] W. Ballik and K. Lake, "The volume of static black holes" >rXiv:1005.1116 

[24] This is a package which runs within Maple. It is entirely distinct from packages distributed with Maple and must be 
obtained independently. The GRTensorll software and documentation is distributed freely on the World- Wide- Web from 
the address http: / /grtensor. org^ 



15 



Appendix A: Ingoing null geodesies 

The arguments presented her^ do not rely on a specific form for H. We are concerned only with ingoing radial 
null geodesies. It is convenient to write the associated null geodesic equation in the form 



dH / -ff(l - f ) 2m\ fdt 



1. r]i 



2 



(Al) 



Along rji, / < and defining 



we have 



Now along Tji, R is strictly increasing to i?+ > 3m and so h{R) eventually becomes negative. Define some fiducial Rq 
so that 

h{Ro) = -a (A4) 

where < a < l/12m. Now (|A4p gives 



i?± ^ 1 ± VI - 12am ^^^^ 
2a 

and since we can always choose a sufficiently small (but not zero) so that Rq < < Rq we eventually have 

dH f dt^"^ 

along Tji for A > some A*. Now write 



so that from (|Aip we have 



Integrating (|A8p we have 



Finally, integration of (|A9[) gives 



T < TTTTT- (A9) 



a(A-A,)-l/r, 



t-i, < -ln(r,a(A- AO + 1). (AlO) 



From (|A10p it follows that t oo for finite A and so T]i is geodesically complete. 



Some of the conclusions in this Appendix (and we expect by now all) have also been obtained by Brien Nolan. 
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2. {b,i+) 

We consider the ingoing geodesies which reach V. within the range {b,i^). Sufficiently close to H, R is strictly 
decreasing and / > with equality holding only on H. We now have 



where we continue to use (|A2[) . Now since R approaches R- < 3m, h{R) is strictly positive and increasing. Again, 
for some fiducial Rq, but now > write /i(-Ro) — a > 0. We now have 

eft f dt^^ 



for A > some A». Now integrating in parallel to the details given in the above case we arrive at 

t-i, > -In ( -— — ^— — — ) . (A13) 



a \T^,a{\^ — A) + 1 

From (lAlSp it follows that i — ^ oo for finite A and so ingoing geodesies that reach T-L in the range (&, i+) are geodesically 
incomplete. This observation is the principle contribution given in KKM. 

3. {-z\b] 

Finally, we consider the ingoing geodesies which reach T-L within the range This includes T]2- Now R is 

strictly increasing and / < with equality holding only on %. We now have 

d'^t , , , f dtV , . , 

d^^'^^^dx) (^^4) 

where we continue to use (jA2p . Now since R approaches i?_ < 3m, h{R) is strictly positive and decreasing. Again, 
for some fiducial Rq < R-, write h{Ro) = a > 0. We now have 

dH f dt'^'^ 



for A > some A». Now integrating in parallel to the details given in the above cases we arrive at 

t-i, < iln ( -— — ^— — — ) . (A16) 



a \r*Q;(A* — A) + 1 

From (jA16p it follows that i — cx) for finite A and so ingoing geodesies that reach T-L in the range are also 

geodesically incomplete. 

Appendix B: Bifurcation two - spheres 

A bifurcation two - sphere is usually discussed in terms of a vanishing time-translational Killing vector (e.g. ilftj). 
Here we define a bifurcation two - sphere in terms of the simultaneous vanishing of both the ingoing and outgoing 
radial null geodesic expansions. First, for clarity, let us review the situation in the Schwarzschild vacuum. As shown 
in J the Kruskal - Szekeres metric can be given as 

ds^ = {2Mfd~s^ (Bl) 

with 

^ T, — ,^r dudv + (1 + Cfd^ll (B2) 
(1 + Lje^^'- 
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where 

C^C{-^) (B3) 

and C is the Lambert W function Trajectories with tangents fc" = e^(l + (constant u = uo,0 and </>) are 

radial nuU geodesies given by 

v{X) = Xe~^ (B4) 
where A is an afhne parameter and we note the expansion 

Trajectories with tangents I" = e^(l + C)S!j^ (constant v — vo,9 and (j)) are radial null geodesies given by 

u{X) = Ae"^ (B6) 

and we now note the expansion 

V„r = — (B7) 

e{l + C) ^ ' 

On the horizons u = and v = Q then v and u are afEne parameters. The bifurcation two - sphere is given by 
u = V = Q and it is uniquely characterized by Vafc" = VqI" = 0. 

The McVittie solution under consideration in this paper reduces to the Schwarzschild - de Sitter spacetime on 
% and so we need a more general construction. We follow jl^ . It is shown there that for all static metrics 

ds-^ = -fde + ^+r^dnl, (B8) 

where / is a polynomial with simple root(s), 

/(r) = (r - a)h{r) (B9) 
where h{a) ^ 0, one can construct regular extensions about r = a via the transformations 

uv = ±{r - a) exp( / 2K'^dr + £), (BIO) 

where the sign depends on how we choose to orientate the u — v axis, £ is a constant and 

= exp{2Kt) (Bll) 



where k is the surface gravity given by 



_ Idf 
2dr 



^0. (B12) 



Note that according to (|B10[) r — r{uv) and r |a 7^ 0, = d/duv. Under these transformations the Killing vector 
7]°' = 5f becomes 77" = {u, ~v, 0, 0) and one recovers the usual definition of the bifurcation two - sphere at u = ?; = 0. 
Note that the specified construction can always be done. However, about a distinct root, say r = 6 7^ a, a new chart 
must be constructed about r — b. 

To calculate null geodesic expansions note that the metric takes the form 

ds^ = K{r)dudv + r'^dnl (B13) 

where 

i^(r).±^exp(-2«/M^d.). (B14) 
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The integration constant has been absorbed into the factor a and again the choice of sign determines the orientation 
of the u — V axis. Trajectories with tangents fc" — -^jpj (constant u — uo,0 and 0) are radial nuh geodesies with 
expansions 

K[r)r 

and trajectories with tangents P = T^fry (constant v ~ vo,0 and (f>) are radial null geodesies with expansions 

K[r)r 

The bifurcation two - sphere associated with any non - degenerate horizon at r = a is therefore characterized by the 
simultaneous vanishing of both the ingoing and outgoing radial null geodesic expansions; VqA:" — Val"" = 0. 



Appendix C: The degenerate case 
Since H < 0, fo > f . In the degenerate case 27m?H^ = 1 and so 

(3m- i?)2(6m + i?) 



/< 



27m?R 



(CI) 



Throughout the associated McVittie solution / < and / = only at the horizon T-L where R ~ 3m. As a result, all 
ingoing and outgoing radial null geodesies have dR/dt > and there can be no bifurcation two - sphere. All ingoing 
geodesies that reach H do so without crossing / = first and as in the non - degenerate case we find that these 
ingoing geodesies are incomplete. The ingoing radial null geodesies in the z — I plane are shown in Figure [T51 
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FIG. 13. As in Figure[5]but for /o = 1/Hf,m — 3\/3 (again we take Hq = 1/3). The limiting geodesies Tji and Tj2 now coincide 
and give Tj. There is no locus f — 0. The dashed curves connecting z = and H and H with z — 1 indicate Sz/dl^ = 0. 

A possible extension of the associated McVittie solution is shown in Figure [14] where we have used the mapping 
function 



i = A(l - a;)tan {^{l - x)'^ 
with A ~ 1.9 and a degenerate Schwarzschild - de Sitter "interior". 



(C2) 
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FIG. 14. As in Figure [121 but now for the degenerate case. 



Appendix D: Ho = 

The case Hq — proceeds in a fundamentally different way than the case Hq > 0. To see this, consider the usual 
background of "dust" so that H = 2/3t. Defining x — R/m and T = t/m the locus f — takes the form 

(Dl) 



and so surfaces of constant t never intersect the locus for t < 2'\/3to, intersect it once for t — 2\/3m and intersect it 
twice for t > 2\/3m. Moreover, every surface of constant R, for 2m < R < oo, crosses the locus once and at finite 
t. This last point shows us that eventually all ingoing radial null geodesies have / > and so eventually dR/dt < 0. 
There is no bifurcation two - sphere in this McVittie spacetime. To integrate the radial null geodesies we continue to 
use (HI]) but replace (gT]) with 

I = i . (D2) 

1 + mH ^ ^ 



The radial ingoing null geodesic equations now take the form 

dz _ (l-z2)2 / 2(1-0 



dl 6(1-/2) \^;(l_^2) 



(D3) 



The integrations are shown in Figurejini Using the same procedures as before, the surfaces of constant R and constant 
t are shown in Figure [ini We find that H is at finite affine distance along ingoing null geodesies but 1+ is at infinite 
afRne distance along outgoing null geodesies (this spacetime is asymptotically flat). 

The conformal diagram and a possible extension onto the Kruskal - Szekeres manifold is shown in Figure 1171 There 
is now a black hole horizon to the future of the singularity and a white hole horizon in the past. These now form a 
wedge, not a straight line. 

Now KKM argue, their considerations motivated by quantum gravity arguments, that the invariant 

A = V,V,i?,;3^5V'V^i?"'^^^ (D4) 

diverges on the horizon. It is difficult to see how this would come about since at the horizon, using conditions ([7]), with 
Hq > 0, the metric tensor along with all first and second order (partial derivatives) are continuous. Using GRTensor 
II [13], and assuming _ff is e C'^, we find 

1440m2 

^hm A = — ^j^(7f2i?3(ii7j2^3 _ + 59^) + UR^ _ QQRm + 65m^) (D5) 
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FIG. 15. As in Figure [5] but for Ho = 0. Note that H is now given by {R = 2m, t = oo). 

,+ 




P R=lm I 

FIG. 16. Gonformal representation of surfaces of constant R (dashed) and constant t (dots) for the case Ho — 0. The locus 
f — from p to i'^ is shown (solid). 



where evaluation along 7y_ is understood. We see no divergence for R > 0. 




FIG. 17. An extension of the McVittie spacetime for Ho = which is nuU geodesically complete. Note that all trajectories 
terminate on the spacelike singularity R = 2m. The loci / = are shown. 



